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The Berger model of perturbative fragmentation of quarks to pions ^ is improved by providing 
an absolute normalization and keeping all terms in a (1 — z) expansion, which makes the calculation 
valid at all values of fractional pion momentum z. We also replace the nonrelativistic wave function 
of a loosely bound pion by the more realistic procedure of projecting to the light-cone pion wave 
function, which in turn is taken from well known models. The full calculation does not confirm the 
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behavior of the fragmentation function (FF) predicted in [1| for z > 0.5, and only works 



at very large z > 0.95, where it is in reasonable agreement with phenomenological FFs. Otherwise, 
we observe quite a different z-dependence which grossly underestimates data at smaller z. The 
disagreement is reduced after the addition of pions from decays of light vector mesons, but still 
remains considerable. The process dependent higher twist terms are also calculated exactly and 
found to be important at large z and/or pt- 
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PACS numbers: 12.38.-t, 12.38.Bx, 12.39.-x, 13.66.Bc 



I. INTRODUCTION 

The fragmentation of colored partons, quarks and glu- 
ons, into colorless hadrons is an essential ingredient of 
any semi-inclusive hadronic reaction, since confinement 
does not allow propagation of free color charges. For this 
reason hadronization is usually considered to be related 
necessarily to confinement specific to the string model Q ■ 
Indeed, the string model of hadron production is rather 
successful in describing data. 

In a typical event of quark fragmentation the mean pro- 
duction time tp of a pre-hadron (i.e. a colorless cluster de- 
veloping afterwards a corresponding wave function) lin- 
early rises with its energy, and the most energetic hadron 
in such event takes about half of the initial quark energy. 
In some rare events, however, the leading hadron may 
take the main fraction z — > 1 of the initial quark energy. 
This process cannot last long, since the leading quark is 



constantly losing momentum, dpq/dt 



where k is 



the string tension. Therefore the production time should 
shrink at z — > 1 as [3|], 



(1) 



Notice that the end-point behavior of the production 
time, tp (X (1 — z), is not specific for the string model, 
but is a result of energy conservation. 

The shortness of the production time is an indication 
that a nonperturbative approach for the production of 
hadrons with large z — > 1 is not really required. In- 
deed, according to ([1]), in this region the hadronization 
time shrinks, i.e. the quark directly radiates a hadron, 
q — > h+q. Furthermore, since the invariant mass squared 
of the final state is M^^ = m\/z+m?q/ (l-z)+py/z(l-z), 
where px is the transverse hadron momentum, at z — > 1 
the initial quark is far off mass shell, and this process can 



be treated perturbatively. This observation motivates a 
perturbative QCD calculations for leading pion produc- 
tion q irq, within the model proposed by Berger 
as is illustrated in Fig. [1] for II annihilation. He found 




FIG. 1: II annihilation with production of two qq pairs. The 
large blob contains gluon radiation by either q\ or q2. Four- 
momenta of particles are shown in parentheses. 

that the fragmentation function of a quark to a pion 
vanishes as (1 — z)^ at z 1, and falls as function of 
transverse pion momentum as Besides, a nonfac- 

torizable, scaling violating term was found to dominate at 
z 1. The shape of z-dependence calculated by Berger 
1] was found to agree well with data after the inclusion 
of gluon radiation cf . Ref. [3| . 

Unfortunately, the calculation performed in [l| missed 
the absolute normalization of the cross section, which 
makes it difficult to compare with data. Moreover, it was 
done in lowest order in (1 — z), therefore it is not clear in 
which interval of z the model is realistic. And last, but 
not least, the calculations were based on the nonrelativis- 
tic approximation for the pion structure function, assum- 
ing equal sharing of longitudinal and transverse momenta 
by the quark and antiquark in the pion. However, the 
dominant configuration of the qq pair projected to the 
pion is asymmetric, with the projectile quark carrying 
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the main fraction of the momentum. 

Here we perform calculations first in the Berger ap- 
proximation, but retaining the absolute normalizations 
and higher powers of {1 — z) (Sect. 3). Then, in Sect. 4 we 
give up the nonrelativistic approximation and project the 
amplitude of qq production onto the light-cone (LC) wave 
function of the pion. For this wave function we consider 
three different models and find reasonable agreement 
with phcnomenological fragmentation functions (FF), 
but only at large z > 0.95. To improve agreement at 
smaller z we add pions originating from decays of p and 
uj mesons, which are produced by the same mechanism, 
and which is depicted in Fig.[T] In Sect. 5 we study higher 
twist contributions, which gives a sizeable contribution 
in semi-inclusive pion production in DIS at moderately 
large and large z. 



II. LEADING HADRONS IN BORN 
APPROXIMATION 

The amplitude of the process II qi + q2 + G 
qi + q2 + qs + <Z4, depicted in Fig. [TJ in the lowest order 
of pQCD is given by, 

A{ll ^ qiq2q3qi) = ^ j/ZH^i, ^i; ^2, As) 

X jj:!'Hpi,^i,h;P2A2,i2;P3,h,i3;P4,^4,i4)- (2) 

Here fci , Ai and k2 , A2 are 4-momenta and helicities of the 
lepton and antilepton respectively; pi, A; and ii are the 
4-momenta, helicities and color indexes of the quarks qi 
and (73 {I = 1,3) and antiquaries 52 and 94 (/ = 2, 4). The 
4-momentum Q = ki + k2. 

The leptonic and hadronic currents in ([2]) read, 

jl!\ki,Xi;k2,X2) = euiJk2)j^V)^Jki) ; (3) 
'^/I'^Hpi: Ai,«i;p2, A2,?2;P3, A3,i3;p4, A4,«4) 
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a=l 

X T^„{piXl,P2h) ju{P3h,P4>^4) 



(4) 



Here Af ^ = (ps+Pa)^ is the gluon invariant mass squared; 
5g — 47ras; a|J^ are Gell-Mann matrices; 

7)<i.(PlAl,P2A2) = MAsfe) [ta^GCQ -pi)7^ 

+ 1,^G{p2 - Q)7p] '^Ai (pi) , (5) 

where G{q) = (q + mq)/{q^ - m^); q = g^7^; wiq is the 
quark mass; and 

jv{p3>^3,P4^4) = Ux^{p4) Iv WAsCPs) • (6) 



III. BERGER MODEL 

In the Berger model [1] the amplitude A of the reac- 
tion II — > nqiq^ is a result of projection of the amplitude 
Eq. ([2]) on the 5- wave colorless state of the (7293 pair 
having zero total spin. The result of the projection is 
proportional to "^Trir = 0) (r is 3-dimensional) with a 
pre- factor ■\/2/rn^ Q, where m,r is the pion mass. Then 
we get, 



A{ll ^ nq.q,) = ^ J^'U^J^^ 



*.(0) . (7) 



Here 



V'i ,^1 A=±l/2 

X jj:!'Hpi,M,ii;p,\i;p,-X,i;p4:,X4:,i4) , (8) 

and the summations X]?=i ^'^'^ Sa=±i/2 sS^(A) 
perform projections to colorless and spinless states of the 
9293 pair, respectively. 

Then we can make use of the relations, 



a— 1 ii 

Y] sgn(A)i;_Ab4)wA(P3) = Ibip + m), 

A=±l/2 IP3=P4=P 

(9) 

and arrive at the following form of the hadronic current. 



j{h} _ 

^ 3%/6Af2 



(il/x + j2/i) , 



(10) 



where 



"A4 {P4) lul5{P + mq)-fyG{Q ~ pi)jf, UAi (pi 

u\Ap4) 75 f7p - —j^T^j (pi) ; 



(11) 



j2p = "A4(P4) 71-75 (P + mq)-f^,G{p2 - Q)-fy WAi (Pl) 

/I r 

'«A4(P4)75 (PlP + "l^)7p 



- {Pit,+Pf,)p-mgpj^+mgQ^, v\Api)- (12) 

Here we applied the algebra of 7-matrices, the Dirac 
equation and 4-momentum conservation, Q — pi+2p+p4^. 
The invariant gluon mass M was defined in (j4]). 

It is convenient to choose the 2-axis along the momen- 
tum pi in the collision cm. frame, and to switch from 
Lorentz 4- vectors (e.g. p^i, p^4, Qfi, etc.) 

to light-cone vectors, {a-^.,a-,a±), where a± = qq ± a^. 
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Since Q — 0, i.e. = Q_, the condition of gauge invari- 



ance, Q^J^p = Q^^J}!''' = 0, takes the form, / 
and J 



(0 



-J 



(0 



ih) 



jih) ^ q^j^gjj ^]-^g product of the lepton and 



hadronic currents can be presented as, 



t(0 j{h) 



(13) 



The typical values of transverse components are \p± \ 



Jl , so 



J 



ih) 



Q 



(14) 



Therefore, the first term in can be safely neglected. 
Then we get. 



T(h) 



3\/6A/2 

2m, 



AP 



q ^ 
Pll 



where ^ = (2 + z)/{2 - z), and 

Q+ Q+ 



(15) 



(16) 



the pion form factor reads, F-^{q^) = (1 + (7^/16k2) 
Then = 2^? . 

These two examples demonstrate the high sensitivity 
of the wave function at the origin to the choice of r- 
dependence. One finds |*.(0)|Lp/l«'.(0)|^,„,, = VS^ « 
5. Therefore, it is difficult to conclude whether the Berger 
model agrees or not with data. 

Another, more realistic option would rely on the pole 
form of the pion form factor, F.^{q^) — k1/{k^ + g-^), 
where Kg — l/6(r^^). Then, 



' chl 

2 1 

l*7r(?^| =-exp(-K3r) 
r 



(21) 



In this case, however, the wave function at the origin is 
divergent. 

The Berger approximation, assuming that the pion 
production amplitude is proportional to the amplitude of 
qq production with equal momenta, would be justified if 
the pion was a nonrelativistic, loosely bound system, i.e. 
TOtt ~ 2mq, 2mq — <C ruq. However, the mean charge 
radius squared is much smaller than the value given by 
such a nonrelativistic model, (r^^) = (4to^ — mj)^^. 

On the other hand, a description of the pion as a rel- 
ativistic bound system has been a challenge so far. 



is the fractional pion momentum. In this approximation 
the invariant giuon mass reads. 



m2(l-z/2)2+p2 
z{l-z) 



= {p+P4y 



Notice that although the second term in pS]) is pro- 
portional to the quark mass (which was assumed in Jj to 
be zero), it should not be neglected. Indeed, after inte- 
gration over pt the interference of the two terms in (fT5)) 
is of the same order as the first term squared. 

In this approximation the fragmentation function gets 
the form, 



64a2 



27mTrmq 



1*^(0)1 



e' + 2(e + l) 



2 zjl-z)' 
{2-zY 

^_ 16 z2(l-z) 
^ "3 



2-z 



{2-zY 



(18) 



The pion wave function at the origin correlates with 
the shape of the parametrization for '^■^{r). In the case 
of a Gaussian parametrization. 



I gauss 



T-3/2 



i-exp(-K?rV2), 



(19) 



the pion form factor has the form, ^^(5^) — 
exp(-gVl6K?). So = 3/8(r2j. 

With a bit more realistic exponential shape, 

l*-WlL-4cxp(-2K2r), (20) 



IV. PROJECTION TO THE LC WAVE 
FUNCTION 

A. Direct pions 

In the light-cone (LC) representation the pion wave 
function depends on the fractional LC momenta of the 
quark, a — p2+/p7r+, and antiquark, 1 — a = p3+/p^+, 
and the relative transverse momentum, k± — apz± — (1 — 
a)p2A^- In this representation the amplitudes, Eqs. ([7]) 
and ([2]), are related as. 



A = 



1 



da 



(27r)3 y ^2a(l - a) 



(22) 

where the qq Fock component of the pion LC wave func- 
tion is normalized to unity, 



da 



d'^k. 



(23) 



In this case the projection of the distribution ampli- 
tude of q2 and qs on the pion LC wave function is more 
complicated that in Berger model {a — 1/2), however it 
can be grossly simplified if one neglects small terms of 
the order of m and k± in comparison with large p2+ and 
P3+ order terms. Then the combination in Eq. ([9]) gets 
the simple form, 

^ sgn{X) v-x{P3)ux{p2) =75 P7T + 0{m,k±) . (24) 

A=±l/2 
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Furthermore, neglecting small terms we arrive at a new 
relation for the hadronic current of Eq. (fTS]) , 



J 



(h) _ 8ffge^^ Sa4ai 1 + (1 - a)z) 



3^6 M2 



1 



w(p4)757±''^(Pi)-(25) 



When the momentum fractions of the quark and anti- 
quark in the pion wave function are a and 1 — a, then 
the invariant mass squared reads, 



m2(l - az)2 + [(1 - a)p^±_ - (1 - z)ki_]'^ 
z(l - z)(l - a) 



(26) 



The light-cone pion wave function can be parametrized 



^^{a,r) = (/)(Q;)?/'(r, a) 



(27) 



If the wave function in momentum representation has a 



monopole form, 'I'^(a, k) oc [k'^ /a{l — a) 



tplr^a) — N Ko{Kry^a{l — a)), 



then 



(28) 



where Kq is the modified Bessel function. Since the mo- 
mentum dependence of ^'^^(a, /c) is poorly known, we 
also performed calculations with a dipole dependent wave 
function in the Appendix, since comparison of the results 
shows the scale of the theoretical uncertainty. 
The parameter k is fixed by the condition. 



dF^q) 



1 



' ch/ 



1.83GeV-^ 



q'=Q 6 

where the pion form factor reads, 

1 





(29) 



(30) 



Thus, the parameter k as well as the normalization 
constant N in ((27)) depend on the choice of function (j){a). 
We consider two popular models (compare with Q): 

Model 1: Standard (asymptotic) shape 0, 



(f>i{a) 



a(l — a) ; 



4 



1^2 



chl 



Model 2: Chernyak-Zhitnitsky model 0, 
(/)2(a) = (/.i(a)(l-2a)2; 



70k^ 

TT 

6 



(31) 



(32) 



(33) 



(34) 



To be specific we will calculate D^* (p^ , z) which is 
the FF of a u quark into -k^ . For the transverse momen- 
tum dependent fragmentation function we have for each 



of these versions (taking into account the longitudinal 
current contribution) , 



dDZ {z,p^^) 



dpi 



= 2 



(1 - zf F^{z,pt) + ez'-^G^{z,pT) 



Here i ^ 1, 2; Ci = 1; C2 = 35/3; 



(35) 



Ft{z,pT) 



da 



(1 - Q:)(/)j(a) l-|-(l-a)z 



X In 



(36) 



Gi{z,pT) = da 



(l-a)20,(a) 



(1 - az)y^af - bi 



X In 



Pt{1 - + m^(l - az)^ 

Kfa{l — a){l — z)^ ; 

Am^qufil ~ azf{l - zfa{l - a) . 



(37) 

(38) 
(39) 



In fact, only the first leading twist term in square 
brackets in ([35)) corresponds to the factorized FF. The 
second term is a higher twist term, whose value (factor 
e) is process dependent, and which is discussed in more 
detail in Sect. IV) below . 

The results of the numerical calculations of the p- 
integrated FF, for each of the three models, are plotted 
as functions of z in Fig. [2l The QCD coupling was fixed 
at as = 0.4. 

The calculated fragmentation functions fall off with z 
only at very large z — > 1, otherwise are rather flat, or even 
rise at small values of z. Such a behavior does not comply 
with data which suggest FF monotonically falling with z 
(lol |. Apparently, the present calculations are missing 
some mechanisms contributing at small z. 



B. Vector meson decays 

One of the processes contributing to the pion spectrum 
should be the production, by the same mechanism shown 
in Fig. [H of heavier mesons which decay to pions. One 
of the most important corrections should come from p- 
meson production, which gives the following contribution 



ADS/-^(z) = 



1 



dz' 



Di iz') + DP, {z') 



(40) 
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FIG. 2: The fragmentation function, Eq. (|35|l . integrated over 
transverse momentum. Solid and dashed curves correspond 
to the models 1 and 2 for the pion LC wave function (see text) 
respectively. 



The bottom integration limit reads, 



2z 



4m: 



(41) 



We assume that Dl^ (z) = 3I?J (z), since p has spin 
1, and that Dp° (z) = ^Dp" (z). 

The w-meson production may also be important. Pions 
from Lo decays should be even softer because of the three- 
particle phase space. The corresponding correction to the 
pion spectrum can be calculated as follows. 



zi 



Z2 



min < 1, 



min < 1, z 



2m^ z 



2m.'?, 



(45) 



We assume that £'"(z) = (z), since the factor of 
3 coming from spin enhancement is compensated by an 
isospin suppression. 

Fig. [3] shows our results for (z) (dashed-dotted), 
(z) and ADu^^ (z) (dotted), and their sum 
(solid). We also plotted the phenomenological D"^ (z) 
(dashed) obtained from a global fit to data As an- 
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FIG. 3: Comparison of the Model 1 (asymptotic shape of the 
pion wave function) with data. The curves from bottom to 
top are: dotted: pions from uj and p decays; dot-dashed: direct 
fragmentation to pions; solid: sum of the three previous con- 
tributions; dashed: phenomenological FF for charged pions 
[Tol l fitted to data at scale p' = 0.5 GeV^ 



dM2^ g(Af2.) /(z, M2^) 



dM2^g{M2-,) 



where 



g{M2.) - ^{Ml,-Aml){n^~Amlml), 



and 



/(z,Af2.)= / '^D-{z')., 



ticipated, the production of p contributes to the softer 
, (42) part of the pion momentum distribution, and does not 
affect its hard part. 

Other meson decays should pull the medium-z part of 
D'!^ (z) further up, but accurate calculation of all those 
contributions is still a challenge. 

Notice that our results have no evolution, since 
the calculations are done in Born approximation. Mod- 
(43) 

ification of the z-dependence by gluon radiation makes 
it softer, closer to data, generating also a evolution. 
These corrections were studied within the Fock state rep- 
resentation in 

The transverse momentum distribution of pions is 
(44) given by Eq. ([55]) . One cannot compare with data the 
mean value of {p^) since it is poorly defined. Indeed, 
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so {p\) is divergent and 



Fi ^ In (pt)/Pt liigl^ PT , so 
depends on the upper cutoff. 

Instead, one should compare with data the pr depen- 
dence. Our results for the pr-distribution of the FF, 
Eq. (pS)) . is depicted in Fig. U] for several values of z. 
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FIG. 4: The transverse momentum dependent FF, 
dDg{z)/dpT, calculated with Eq. (|35j and the Model 1 for the 
production of direct pions. Solid, dashed and dotted curves 
are calculated at z = 0.75, 0.85 and 0.95 respectively. Data 
from df] at IF^ > 350 GeV^ are renormalized for a better 
comparison with our results. 

It might be too early to compare these results with 
data, since we did not include yet the gluon radiation, 
intrinsic motion of quarks in the target, and decays of 
heavier mesons. Nevertheless it is useful to check whether 
the calculated px dependence is in a reasonable accord 
to data. Notice that the data depicted in Fig. 0] are in- 
tegrated over a rather large z-bin, 0.4 < z < 1. The lat- 
ter causes a considerable mismatch in normalization (see 
Fig. [3]), so we renormalized the data [fjj to be able to 
compare the shapes, which then are in reasonable agree- 
ment. 



V. HIGHER TWIST TERMS 

The last term, in square brackets in Eq. psp . is a 
higher twist effect. It does not vanish at z ^ f, but 
is suppressed by powers of Q. We neglected corrections 
of the order of {p"^) / {zQ'^) , which are important only at 
small z. 

This higher twist term breaks down the universality of 
the fragmentation function, since the factor e depends on 
the process. For e~^e~ annihilation it is given by. 



I + cos^0 



(46) 



where 9 is the angle between the direction of II collision 
and momentum pi in the cm. frame. 
For deep-inelastic scattering it reads. 



1-y 



2(1-2/) +2/2 



(47) 



where y = is 4- momentum of the virtual pho- 

ton; I is 4-momentum of the initial lepton. 

The relative contribution of the higher twist term is, 



Ri{z,pT) = 4e 



1-zJ Q'^Ff{z,pT) 



(48) 



where subscript i denotes the number of the model used 
for the LC pion wave function, and Gi, Fi are defined in 

While the relative value of the nonfactorizable higher 
twist term is expected to be vanishingly small in II an- 
nihilation, it might be a sizeable effect in SIDIS, usually 
associated with medium to large values of . The rela- 
tive correction, Eq. is plotted in Fig. [5] as function 
of Pt, for = 2.5 GeV^ and several fixed values of z. 
Solid and dashed curve correspond to the models 1 and 




Pt (GeV) 



FIG. 5: The relative higher twist correction to the FF of a 
quark in DIS as function of transverse momentum for fixed 
values of z = 0.75, 0.95, and = 2.5 GeV^ The solid and 
dashed curves come from calculations with Model f , Eq. (|31|l . 
and Model 2, Eq. (|33|l . respectively 

2 for the LC pion wave function, respectively. Although 
the higher twist term is relatively small for forward frag- 
mentation, it becomes a dominant effect at > 1 GeV^. 

The corresponding higher twist correction to the pt- 
integrated FF reads. 



R^{z)=Ae■ 



/ z V s;;°dpiGi{z,pT) 

Q2 j^dplFf(z,pT)' 



(49) 



The factor {p^) is divergent and depends on experimental 
kinematic cuts. Therefore one should rely on its value 
specific for each experiment. 

Apparently, a direct way to see the higher twist con- 
tribution in data is to study the behavior of the FF. 
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However, such data at sufficiently large z are not avail- 
able so far. Therefore, we try to extract the higher twist 
contribution from the z-dependence. To do so we first 
fit data at moderate values z < 0.65 where we do not 
expect a sizeable higher-twist corrections, with the stan- 
dard parametrization D'^{z) = Nz"(l — z)^. We use 
data from the HERMES experiment [12]. We added the 
statistic and systematic errors in quadratures. The data 
are corrected by subtraction of the contribution from 
diffractive vector mesons, j*p 7rp, which is another 
higher twist contribution (see section IVI[) . We found 
a = -1.24 ± 0.04, /3 = 1.5 ± 0.07, N = 0.88 ± 0.07. 
The data divided by this fitted z — dependence are de- 
picted in In Fig. [6] We compare this data with the rela- 
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FIG. 6: Hermes data T^] for multiplicity of charged pions 
produced in DIS on a proton, corrected for decays of vector 
mesons. The data points are divided by the fit to the data 
at z < 0.65 (see text). The curve corresponds to Ri{z) + 
1 calculated with Eq. (gg at = 2.5 GeV^ and (py) = 
0.25 GeV^ 

five contribution of higher twists Ri{z), Eq. calcu- 
lated at — 2.5 GeV^ and with the measured value of 
(p^) fa 0.25 GeV^ [l2|. Our results agree with the data 
reasonable well. 

An attempt to see the higher twist effects in nuclear 
attenuation data was made in p^ . They found higher 
twist corrections of similar magnitude. 

Notice that other sources of pions, like decays of heav- 
ier mesons produced via the same mechanism, are impor- 
tant for leading twist part. However, they also supply the 
cross section with higher-twist terms. Nevertheless, we 
assume that these corrections affect the ratio much less 
than the cross section. 



VI. HINTS FROM TRIPLE-REGGE 
PHENOMENOLOGY 

The factorized part, Eq. of the cross section of 

pion production in II annihilation, is the same as in deep- 
inelastic scattering (DIS), where it can be compared with 
the expectations of the triple-Regge description, illus- 
trated in Fig. [71 The inclusive cross section at fixed z 
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FIG. 7: Virtual photoproduction of a pion via Reggeon ex- 
change. The projectile quark from the photon fluctuation 
picks up an antiquark, produced either from the vacuum or 
perturbatively (see Fig. [T} , and they form a pion. 



is energy independent (Feynman scaling), and at hxed 
energy and 1 — z <C 1 depends on z as. 



da{'j*p ttX) 
dz dp^ 



(x(l-z)". 



(50) 



where z equals to Feynman in the triple-Regge kine- 
matic region. 



Xf = \ l ) (1 - XBj) 



(51) 



and XBj is the Bjorken variable. 

The exponent in ([50]) is related to the parameters of 
the Regge trajectories involved, 



n = 1 - 2a]B,{pj 



(52) 



Here a][i{p\) is the trajectory of Reggeon M. The rapid- 
ity interval. Ay « — ln(l — z), covered by the Reggeon 
is not large for the values oi z ^ 0.9 under discus- 
sion. Therefore the pion Regge pole should dominate, 
since it has large coupling to nucleons. In this case, 
(^T^iPr) ~ "Qf^rPT; where a'^ ~ \ GeV^^. Thus, 



1^ = 2a'Apl) -1-5. 



(53) 



Here we rely on the value (p^) « 0.25 GeV^ measured 
in both HERMES [12] and EMC Till experiments. The 
value of the exponent given in Eq. (j53p agrees quite well 
with data. Although our calculation confirmed the value 
n = 2 found in 'l'], the inclusion of gluon radiation re- 
duces the exponent n down to the value observed in data 
i. 

Notice that the z-dependence presented in Eqs. ([50)1 - 
([5^ changes at very small 1 — z <C 1, and becomes rather 
flat. Indeed, we assumed that the invariant mass squared 
of the excitation X is sufficiently large, s{l — z) ^ for 
the Pomeron to dominate in the bottom leg of the triple 
Regge graph in Fig. [71 However, this condition breaks 
down at very small 1 — z and Reggeons with aK(0) = 1/2 
dominate in the bottom leg. Another assumption we have 
made, pion dominance in the i-channel exchange, is also 
violated when the rapidity interval ln(l — z) becomes very 
large. Then Reggeons with a higher intercept Q!iR(0) = 
1/2 become the dominant contribution. Thus, the end- 
point behavior has the same power dependence, Eq. (IST 
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but with a different exponent, 



VII. SUMMARY 



q;r(0) - 2a][i{p^) 



(54) 

Tfius we arrive at the remarkable conclusion that the 
FF, which falls steeply with z, levels off at very small 
1 — z ^ 1. This behavior, dictated by the triple- Regge 
formalism, is more general than perturbative calcula- 
tions. One may wonder why this end-point feature is 
absent in our calculations. What has been missed? No- 
tice that we did not care about the fate of the recoil 
quark (74 in Fig. [H which was justified by the condition 
of completeness. However, if the target excitation X has 
a small invariant mass, it affects the probabilities of dif- 
ferent final states of 94. 

The triple-Regge approach also indicates as an addi- 
tional source of a higher twist contribution, which is spe- 
cific for semi-inclusive DIS (SIDIS), the diffractive inclu- 
sive process 7*p pX. The pr-integrated cross section 
corresponding to the triple-Pomeron graph can be pre- 
sented in the form, 



da{j*p pX) 
dz 



p 



167r 



(l-z)|ln(l-z)| {aZ? 
da{-i*p ~> pp) 



dp^p 



(55) 



where G^p{Q) = 3.2mb/GeV^ is the effective triple- 
Pomeron coupling, extracted from the fit [lij to data on 
pp pX. Here we neglected the transverse size of the 
qq dipole projected to p, since it is small, and the 

Pt dependence of the bare triple Pomcron vertex, since 
it is very weak [l^ . All the cross sections in (|55)) should 
be taken at a cm. energy squared s' = so/(l — z), where 
So = IGeV^. 

The z-distribution of the produced p*^-mesons strongly 
peaks at z ^ 1 (as any diffractive process should) and 
their decays feed the effective FF D'^{z), 



We performed calculations for the Berger perturbative 
mechanism T] of quark fragmentation into leading pions, 
keeping all the sub-leading terms in powers of (1 — z) 
and all the coefficients. Our results can be summarized 
as follows. 

• We performed a full calculation of the quark FF 
including higher twist terms within the Berger ap- 
proximation. However, we concluded that the ap- 
proximation of a nonrclativistic pion wave function 
is unrealistic and brings too much uncertainty to 
the results of the calculation. 

• We projected the produced qq pair distribution am- 
plitude to the light-cone pion wave function. For 
the latter we employed two popular models: (i) 
the standard asymptotic shape PT|) : (ii) Model of 
Chcrnyak-Zhitnitsky p3|) . Both models lead to a 
z-dependence quite different from the one inferred 
from data. Only at z > 0.95 our calculations agree 
reasonably with data (both the shape and value), 
but greatly underestimate data at smaller values of 



• Remarkably, the main amount of pions produced 
in quark fragmentation are not produced directly, 
except the most energetic ones with z > 0.95. This 
fact should be taken into account in models em- 
ploying perturbative hadronization 

• Searching for ways of improving the description of 
data we added pions originated from decay of light 
vector mesons p and lo. Although this contribu- 
tion pulled up the production of pions at medium 
to large z, apparently some contributions are still 
missing. That may be production and decays of 
heavier mesons, which are difhcult to evaluate. 



diff 



1'P 
"tot 



dz' da{-f*p p°X) 



(56) 

Here ^ and z,„i„ are defined in (|^T|l . Due to color trans- 
parency the amplitude of rho production is inversely pro- 
portional to Q^, therefore a{'-f*p P^X) oc l/Q**. On 
the other hand, the total virtual photoabsorption cross 
section is a^^f cx l/Q^ (Bjorken scaling). Therefore, the 
diffractive contribution to the effective FF g ^ tt is a 

higher twist effect, IaD?,^''* (z)] cx l/g^. 

The elastic production of vector mesons, 7*p — > Vp 
certainly also contributes to inclusive pion production, 
and is also a higher twist effect. It can be evaluated 
using Eq. and a delta function for the z'-distribution 
of pro duced vector mesons. However, in some cases, like 
in [13, this contribution has been removed from data. 



• We also performed a full calculation for the higher 
twist term originated from the longitudinal current 
contribution. It overcomes the leading twist term 
at large z and/or large transverse momenta. 

• A new higher twist contribution to pion produc- 
tion is found. It is related to decays of diffractively 
produced vector mesons. 

It worth reminding that our results for the FF at large 
z > 0.9 should be compared with a phenomenological 
one with precaution. First of all, data at such large z are 
scarce and different parametrizations [l^, [H, [13 differ 
from each other considerably. Second of all, our FF is 
calculated in the Born approximation. Evolution (gluon 
radiation) may considerably change the shape of the z- 
dependence 
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Appendix A. 



N 



10 



10 



DIPOLE FORM OF THE PION LC 
WAVE FUNCTION 




verse momentum dependent part of the LC wave function 



"^Tr{a, k) cx 



a(l — a) 



In impact parameter rep- 



resentation it takes the form (compare with ([27])), 

-^^{a^r) = N (j){a)^ya{l - a) rKi{Kry/a{l - a)), 

(A.l) 

In this case we can still employ Eq. (|35p for the frag- 
mentation function, but with a new form of function 



Fi{z,p) ^ J da 



(1 — a)(f)i{a) 1 + (1 — a)z 



\ ~ az 



d.,{a-2ei) I ai + ^/ af - bi 
Ui — 2di H , _ m ' 



where di = Kfa{l — a){l — z)^; — m^{l — azY 
Parameters Ci and Ki in (j35[) also get new values, 
Model 1: asymptotic shape. 



(A.2) 



Nl = 



9k? 

27r ' 
36 



Ci = 3. 

Model 2: Chernyak-Zhitnitsky shape, 



(A.3) 



FIG. 8: Fragmentation functions for direct pions calculated 
with pole, Eq. (^7)) (solid curves), and dipole, Eq. (|A.ip 
(dashed curves), parametrization for the transverse momen- 
tum dependent part of the LC pion wave function. Labels 1 
and 2 indicate the model used for the longitudinal momentum 
dependence of the pion wave function. 



. ^ 105^^ 
^ 27r ^ 

108 



Co = 35 



(A.4) 



To see the sensitivity to the form r-dependence of 
the LC wave function of the pion we also performed 
calculations with the dipole parametrization of trans- 



The results of numerical calculations are depicted in 
Fig. [8] in comparison with calculations performed with 
the pole parametrization for the pion wave function. 
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